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Abstract 

We calculate the cross section for the production of pairs of scalar lepto- 
quarks (sleptoquarks) in a supersymmetric Eq model, at hadron colliders. 
We estimate higher order corrections by including vr^ terms induced by soft- 
gluon corrections. Discovery bounds on the sleptoquark mass are estimated 
at collider energies of 1.8, 2, and 4 TeV (Tevatron), and 16 TeV (LHC). 

PACS numbers: 14.80.-j, 11.25.Mj. 
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A natural explanation for the proliferation of fermions, and their pattern of masses 
and mixing angles, is to assume that the current elementary particles are composite. In 
various "preonic" models quarks and leptons have some common constituents. Among 
the particles contained in these models, leptoquarks are exotic particles which have both a 
nonzero leptonic and baryonic number. They also appear in many extensions of the Standard 
Model 0. They can decay directly into a quark-lepton pair, which is a new feature, since 
there are no quark-lepton-boson interactions in the Standard Model. Here we are concerned 
with the leptoquarks that come from a supersymmetric grand unified Eq theory (the low 
energy limit of an ® Eg heterotic string theory [^), where they — together with their 
supersymmetric scalar partners, sleptoquarks — are fractionally charged color triplets. For 
constraints on the parameters of these superstring-inspired models, see, for example, Ref. 
^ or for a more extensive review, Ref. . 

In principle, the best place to look for a leptoquark signal would be at ep colliders 
since they can be produced directly via the lepton-quark-sleptoquark coupling (called from 
now on the Yukawa coupling). However, for small couplings (or large leptoquark masses), 
hadron colliders seem more appropriate since leptoquarks can still be pair produced through 
their strong interactions for arbitrarily low Yukawa couplings. Experiments at LEP have 
already imposed a lower bound of 45 GeV on the leptoquark mass [0 . Results from HERA 
have ruled out masses below 180 GeV for a Yukawa coupling of electromagnetic strength, 
whereas searches at pp colliders have set these bounds to 113 GeV (for a branching ratio 
BR= 1), or 80 GeV (BR= 0.5) 0. More recently, experiments both at HERA and the 
Tevatron have strengthened these bounds reaching with a 95% confidence level , 240 GeV 
for HERA [|1^] and respectively 133 GeV and 120 GeV at the Tevatron []1T|. Note that 



HERA bounds depends on the value of the Yukawa coupling whereas the hadron collider 
results are somewhat insensitive to this parameter. 

Low energy data {e.g. atomic parity violation, . . . ) also impose very strict bounds 
on leptoquark masses. Leurer has updated a previous analysis by Buchmiiller and Wyler, 
and obtained bounds that restrict the leptoquarks (with Yukawa coupling equal to the 
electromagnetic strength) to have masses larger than 600 GeV or 630 GeV for leptoquarks 
that couple to RH quarks, and above 1040 GeV, 440 GeV, and 750 GeV for the SU{2)w 
scalar, doublet and triplet leptoquarks, respectively, that couple to LH quarks [|l^. These 
are "unavoidable" bounds, in the sense that they are independent of the following three 
assumptions, which are used to circumvent other constraints: (1) leptoquarks do not also 
couple to diquarks, (2) leptoquarks couple chirally {i.e. to one quark chirality at a time), and 
(3) leptoquarks couple diagonally {i.e. to a single fermion generation at a time), but they 
remain heavily dependent on the strength of the Yukawa coupling. Finally, some model- 
independent bounds are also discussed by Davidson et al. in Ref. |0. Hence, it is quite 



possible that leptoquarks could escape detection at ep coUiders but could still be detectable 
at hadron colliders. Also, in practice, the energy available in hadron colliders is greater 
than that in ep colliders, which would in principle extend the search to larger values of the 
leptoquark mass. 

The purpose of this brief report is to reanalyze the production of pairs of sleptoquarks 
in hadron colliders |p!3| , p!^ in view of the most recent data and experimental situation. In 
previous work, we calculated the total cross section at center-of-mass energies of 2, 16, 
and 40 TeV (corresponding to the Fermilab Tevatron, the CERN LHC, and the late SSC, 
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respectively), for sleptoquark masses up to 400 GeV. We extend here our calculation of 
the cross section to sleptoquark masses of up to 2600 GeV, and we add the results for 
the (current and newly proposed) Tevatron energies a/s = 1.8 TeV and 4 TeV. (We keep 
the SSC energies in order to compare with our previous results.) The soft-gluon tt^ term 
corrections are included. We also have taken special care to reduce uncertainties in the 
numerical calculations by performing a more reliable numerical integration (using VEGAS 
\T^), and by taking a more recent set of distribution functions. 

Let us briefly review the model and our calculations. (More details are given in Refs. 
n|, and references therein). In the supersymmetric Eg model, the leptoquark D is an 



exotic colored particle which lies in the 27 fermionic multiplet, and the sleptoquark that 
we consider is its scalar superpartner. If we restrict our study to the first generation of 
fermions, then the Yukawa interactions take the form: 



Cy = XlD'* (clUl + i^LdL) + XnDelul + h.c 



(1) 



where c denotes the charge conjugate state, and D is the scalar superpartner of D. We 
assume the A's to be independent and arbitrary. We analyze the pair production of slep- 
toquarks which arise from two subprocesses: (1) quark- ant iquark annihilation {ur + — >■ 
D+D* and ml+m^ D^^+D"), and (2) gluon fusion {g+g D+D* and g+g D^*+D^). 
The subprocess (1) occurs in s-channel (through the exchange of a virtual gluon) and in t- 
channel (virtual electron). The subprocess (2) arises via color gauge interactions from the 
trilinear term gDD in the s-channel (through the exchange of a gluon) and in the t- and 
u-channels (exchange of virtual sleptoquarks), and from the quartic term ggDD in which 
two gluons annihilate to produce directly a pair of sleptoquarks. Note that, to lowest order, 
the gluon fusion process does not depend on the Yukawa interactions. The details of the 



calculation of the amplitudes and the cross sections are given in Refs. JT^ . 

We also estimate vr^ terms, which are the soft-gluon corrections that arise from the 
regularization of either coUinear or infrared singularities, when a timelike momentum transfer 



is involved in the process |T6|. The first-order corrections to a subprocess involving massless 
particles contain an infrared singularity of the form 
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where //^ is the renormalization point. The vr^ term is generated by the last term. 



Re 



(2) 



(3) 



Note that vr^ terms arise only when is timelike. Unless there is a suppression of the 
TT^ terms, their contribution is large and cannot be neglected. It is usually expressed by 
means of the so-called K-factor. In general, we expect vr^ terms to appear to all orders in 
as- Low-order results in many QCD processes indicate that the summation of the ensuing 
large corrections into an exponential could take place, as is the case for electromagnetic 
interactions. It is important to remark that not all of the first-order correction diagrams 
lead to a tt^ term. (The detailed investigation of the soft-gluon corrections leading to a K- 



factor through vr term is described in Ref. [14|, for qq and gg subprocesses. The analogous 



study for other leptoquark processes is given in Ref. |jT6| 
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We now display the contributions to the total cross section for each of the two subpro- 
cesses mentioned previously, including the if-factors coming from soft-gluon corrections. For 
the Born term only, one just has to set = 1. The contribution from the qq annihilation 



is: 



1 

<yAB = I dy 



dxKABaAaBiyAB Gu/p{x, Q'^)Gu/p{y, Q^) + (1 ^ 2,p ^ p) , (4) 



with Xmin = 4(M|,+p|,)/s (we takepT = 10 GeV), and Q"^ = s/2 {s = xys). Here AB stands 
for SS, YY, or SY {i.e. purely strong, purely Yukawa, and mixed subprocesses, respectively) 
where as and ay are the QCD and Yukawa couplings, respectively. We set ay = ctem in the 
numerical calculations. Finally, Ga/h{x,Q'^) stands for the distribution function associated 
to the parton a in the hadron h, with scaling variable x and momentum scale Q^. We 
use the parametrization set B2 of Morfin and Tung which is based on more recent 
nucleon data than the parametrization we used in our previous work. Defining x = ^d/^ 
and t] = \Ji — 4:Ml)/s, the purely strong, mixed, and purely Yukawa contributions are, 
respectively. 
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The /^-factors are found to be 

Kss = Ksy 



1 + 7Tas ( Cp - -Ca 



i-y- 2x 
l + ri-2x, 



1 — -Tras, and Kyy = 1. 
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(7) 



(8) 



Here = | and Ca = 3 are the Casimir operators for the SU{3) fundamental and adjoint 
representations, respectively. (Later on, Kss and Ksy will be denoted generically as Kqq.) 
For the gluon fusion subprocess, the integrated cross section is 



(J. 



99 



with 



cr, 



99 



6s 



dy f dxKggalaggGg/p{x,Q'^)Gg/p{y,Q'^), 



g + ^xjr/ + (4 + x)xln 



(9) 



(10) 



The values for Xmin, Pt, are the same as in 
larger correction to the Born term. 
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1 + ^-n-asCp 



Here, the i^T-factor induces a much 

(11) 



1 + gTras + 
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The cross sections and (|T^) are similar to those in (for qq — > gg, gg — > qq), except 

for a factor. For the Born term, one just replaces if by 1 in (Q) and (|Tl|) . Fig. |l] displays 
the total cross section, that is, the sum of gg fusion and qq contributions, with and without 
corrections. 

Higher order effects (only the tt^ terms are included here) are easy to estimate. From 
(1) the TT^ corrections for the quark- ant iquark annihilation are small and negative [i.e. 
they slightly suppress the cross section). For sleptoquark masses in the range 100 GeV 
< M^i < 2600 GeV, the strong coupling constant as lies between 0.13 and 0.09, which 
corresponds to K^q = 0.932 and 0.953, respectively. Here the cross section is suppressed 
by 7 to 4%, depending on the relative importance of the different channels. For the gluon 
fusion process ([Tl|) , the cross section is significantly increased with Kgg = 1.272 and 1.189 
for as = 0.13 and 0.09, respectively. From Fig. |I], we see that the vr^ term leads to a 
slight increase of the cross section for the small values of the leptoquark mass, and that 
the signal is suppressed for large values. This is due to the fact that, at high energies, the 
gluon fusion subprocess-which undergoes bigger corrections- dominates in the region of low 
leptoquark mass. The value of the mass for which the two curves (Born and corrected) 
intersect increases with the center-of-mass energy. These corrections may be crucial if the 
values of ay and are such that we are close to the detection threshold. It should be 
noted that the i^'-factors found here are only part of the first order corrections. However, 
whereas one expects Kgg to be of the same order of magnitude as the rest of the first order 
corrections, it appears that Kgg is the dominant part of the gluon fusion corrections. 

We expect the final leptoquarks to decay into a quark and a lepton. The most interesting 
signature for DD* pair is the production of 2 jets +/"'"/". The standard background would 
come from pp — QQ, where Q is a heavy quark which decays semileptonically into 2 jets 
+z/z7, and will in general involve some missing pt- It is easy to see that the subprocesses 
involved in QQ production get similar vr^ corrections (i.e. a factor Kgg for qq QQ and 
a factor Kgg for gg —>■ QQ). This could slightly modify the signal-to-background ratio 
in a manner which depends on parameters such as the masses of sleptoquarks and heavy 
quarks and the relative importance of each subprocess. For example, if gluon fusion is 
the dominant process in both sleptoquark and heavy quark pair production, the signal-to- 
background ratio will not be affected by the soft-gluon corrections, but the total cross section 
will be appreciably enhanced. 

Finally, we discuss the detection possibilities, taking into account the expected lumi- 
nosities at the various colliders considered in this work. In all cases, we assume that the 
discovery limit is 20 events, and we base our estimates on the corrected cross sections {i.e. 
including i^-factors) . We have obtained the results of Fig. |T| by putting ay = aem- For the 
sake of comparison we give also an estimate for the case of a very small Yukawa coupling as 
well, by setting ay = 0. Experiments at Fermilab expect to gather about 100 pb~^ during 
the current Tevatron run (for which ^/s = 1.8 TeV), which should last until the middle 
of 1995. With the discovery limit mentioned above, the mass reach would be 184 GeV (if 
ay = rather than aem, then it is 181 GeV). If we assume a luminosity equal to 500 pb~^ 
at y/s = 2 TeV, then this number would increase by an appreciable amount: up to 251 GeV 
(247 GeV for ay = 0). The design luminosity at LHC is expected to be 10^"^ cm~^ s~^. A 
run of one year would then give a luminosity of 100 fb~^, providing a search limit of 2110 
GeV (or 2070 GeV for ay = 0). Fermilab has recently started discussing an upgrade of 
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the center-of-mass energy to 4 TeV, and the numbers used by different people as expected 
luminosities are 1 fb^^ and 10 fb^^. The discovery of 20 events would be reached for masses 
not exceeding 434 (428) GeV and 597 (590) GeV, respectively (numbers in parentheses are 
for ay = 0). At the SSC, with a hypothetical luminosity of 10 fb~^, the same discovery 
limit would lead to a mass reach of 2565 GeV for ay — asm, and 2525 GeV for ay — 0. 
Actually, changing from ay = aem to ay = has an overall effect of only a few percent on 
the discovery limits. This means that the strong processes dominate even at high energies 
where the ratio Q;em(Q^)/Q'5(Q^) is not that small. 

We are indebted to N. Hadley for information on the experiments being carried out at 
the Fermilab Tevatron, P. Labelle for discussions about our numerical calculations, and D. 
London for reading the manuscript. L.M. is grateful for the hospitality of the Service de 
Physique Theorique (CEA-Saclay), where part of this work was done. This research was 
supported by Fondation A.F.D.U. Quebec (G.S.), the Natural Sciences and Engineering 
Research Council of Canada, and by the Fonds pour la Formation de Chercheurs et I'Aide 
a la Recherche du Quebec. 
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FIGURES 

FIG. 1. The total cross section for pair production of scalar leptoquarks at pp colliders with 
a center-of-mass energy of 1.8, 2, 4, 16 and 40 TcV versus the mass of the leptoquark. The Yukawa 
coupling is taken to be ay = Oem- Solid (broken) lines represent the Born (vr^-corrected) cross 
sections. The bottom- left axis are valid for 1.8 and 2 TeV, and the top-right axis, for 4, 16 and 40 
TeV. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9406321vl 



